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INTRODUCTION 


Many scientific investigations have “been pursued in the past few years 
concerning the radiation characteristics of slotted antennas immersed in a 
plasma. This research has, of course, been prompted by attempts to explain or 
predict electromagnetic effects which may occur during the reentry of space 
vehicles. 

The analysis of the problem is conceptually simple, i.e., the wave equa- 
tion is solved, and all unknown coefficients are found by applying the electro- 
magnetic boundary conditions at the air-plasma interface and at the aperture. 
However, unless idealized models of the vehicle are chosen, the evaluation of 
the complicated solutions of the wave equation becomes a major obstacle. In 
order to avoid these problems, models chosen by most investigators have been of 
simple geometry, namely, the coated slot on the flat ground plane and the slot 
on a conducting cylinder. The former is applicable to the reentry problem pro- 
vided the dimensions of the vehicle are much greater than a wavelength, and the 
latter is a reasonable choice if the wavelength is comparable to vehicle 
dimensions. 

Radiation pattern computations constituted the earlier work in this field. 
The slot on the flat ground plane, coated with a gyro-plasma was analyzed by 
Hodara (refs. 1 and 2), and the coated cylinder was originally solved by Wait 
(ref. 3). Others (refs. 4, 6, 7> 8, and 10) have also examined these 

geometries from different viewpoints or to extend the results. The final 
expressions for the patterns reduce to relatively simple forms because of cer- 
tain asymptotic expansions which may be used in the far field. At least for 
plane geometry, patterns may be computed easily without the aid of electronic 
computing machines. The results are valuable because changes in the angular 
distribution of energy can be predicted. 

However, the radiation patterns do not give a complete description of the 
problem. For instance, if the aperture is backed by a waveguide, the amount of 
power reflected back into the guide cannot be established by the patterns alone. 
Such information requires knowledge of the input admittance which follows from 
a study of the near fields. 
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Recently, excellent work on this subject has been done by Villeneuve , . 
(ref. 11 ), Galejs (refs. 12 , 15 , and l4) and Compton (ref. 15 ) • All three con- 
sidered a waveguide opening onto a flat ground plane; however, the differences 
in the formulation and results made each effort a contribution. 

The only published near-field results for the cylinder is Wait's treatment 
of the thin half-wave axial slot radiating into free space (ref. 3, p. 48) . 

|_Jji this paper, the near fields of both geometries are investigated via a 
transfer function which algebraically relates the unknown aperture fields to 
those in free space. This transfer function is a matrix of computable numbers 
which are specified entirely by the plasma parameters, regardless of whether 
the plasma is homogeneous or inhomogeneous. It is important to note that the 
matrix is independent of the aperture excitation. Therefore, when a given 
plasma condition is analyzed, any aperture problem can be solved with no more 
degree of difficulty than the same aperture radiating into free space. The 
other advantage of this type of formulation is that better near-field approxi- 
mations may be made without repeating the plasma solution. 

To illustrate the use of these techniques, the admittance of long slots 
on cylinders and planes is computed. The results are compared, and the effects 
induced by the inhomogeneous plasma are noted. Although no computations for 
finite apertures are given, expressions are derived. 

An appendix is also included which outlines an exact treatment of this 
parallel plate waveguide opening onto a flat ground plane. 


RADIATION FROM SLOT ANTENNAS ON A FLAT GROUND PLANE 


I. Infinite Slot with TE X Polarization 

The model shown in figure 1 consists of a conducting ground plane onto 
which is cut a long radiating slot excited so that H is directed along the 
slot. The fields are assumed to be uniformly distributed along the slot, and 

the exponential function e*^ 0 ^ is chosen to describe the time dependence. The 
complex permittivity of the plasma is assumed to vary only in a direction 
normal to the ground plane. 

The scalar wave equation which describes the magnetic field in the inhomo- 
geneous plasma is 


^ 2 H x I (y^ z ) . ^ 2 H x I (y,z) 1 de ^ 2 x I (y, z ) 2 ( If „\ 

— ^ + — £5 ? + <DTi 0 e(z)Hx(y,z) - 0 


( 1 ) 


The above equation is partially solved by expressing as a Fourier inte- 

gral, i.e., 



Hx 1 ^ 2 ) = ^ f 2 x (ky, z) e ^ k y^cLky (2) 

—00 

which leads to the total differential equation. 

- * ft £ ♦ (5) 

The dielectric constant e is an arbitrary function of position; therefore 
equation ( 3 ) must be solved by numerical means. However, since equation (3) is 
a second-order differential equation, an initial value of H x and its first 
derivative must first be specified before any numerical work can be performed. 
These starting conditions are found by constructing functional solutions of Ey 
and H x in regions I and II, and requiring continuity of the components at the 
interface z = z Q . 

In the plasma, the fields are of the form, 

- s / S x( k y> z ) e31Cyyai y 

^ —oo 

And, in unbounded free space; 
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The preceding choice of roots is required in order to assure outgoing, damped 
waves at z = oo. 

The boundary conditions at z = z Q requires that 


H x I (y^o) = H x n (y,z 0 ) 
E y I (y^ z o) = % n (y^ 0 ) 


( 7 ) 


Therefore, the normalized initial value of H x and its derivative is: 





k 0 < |ky| \ (8) 


ko > |k 0 | 


If the differential equation (3) is now divided by Ep^ky^e^ zZ °, normalized 
values of H x can be numerically computed at z = 0 . 

The normalization constant, Hrp^ky^e^ zZ ° is found by completing the 
boundary conditions at z = 0. 

It is of interest to note the plasma problem has been solved without any 
detailed thoughts given to the aperture or the method of feeding the aperture. 

The transforms of the external field at the aperture are: 


H x (k y ,°) = g(k y ,0)B r (k y ) e j k z Z o 

Mv 0 ) = ;:, r " <3 e (o) 6 , ( k y> 0 )M k y) eJkz2 


k 0 y 0 


J 


(9) 


In matrix form, equation (9) becomes: 



( 10 ) 
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where the matrix elements are known numbers which are defined through equa- 
tions. (9)- Equation (10) is nothing more than a transfer function which 
relates the unknown fields at the aperture to the unknown free— space Fourier 
coefficient Hq>^ky). Therefore when the matrix elements are specified, the 
coated antenna is no more complicated to solve than the uncoated problem. 


If the aperture is fed hy a parallel plate waveguide, excited such that 
only TE X modes exist in the guide, the formal solution follows hy con- 
structing interior (z < 0) fields and requiring continuity of H x and Ey at 

z = 0. The interior solution will contain one unknown, namely, the reflection 
coefficient; therefore, the houndary conditions should, in principle, give hoth 
r and Hr^ky) . Such a general formulation is outlined in appendix I. 

The general equation developed in appendix I has not heen programed; how- 
ever, a first-order estimate of F and Bp can he achieved hy intelligently 
guessing the electric field at the aperture. If a TEM mode is incident on the 
aperture, a reasonable guess for Ey(y,0) is 

Ey(y,0) = ^ (11) 

A construction of the field solutions in the guide to relate the aperture 
voltage V G to the reflection coefficient is the next natural step. Use of 
equation (10) would then give a first-order estimate of hoth F and Bp. How- 
ever, it is more conventional to solve for the external admittance of the guide. 
This can also he done through use of equation (10) . 

The Fourier transform of equation (l4) is 

Ey(ky,o) = V 0 e 72 1 -^ (12) 

k v £■ 

y 2 

Therefore, equation (10) may he expressed in terms of V Q as follows: 



The power per unit length radiated hy the aperture is 

oh/ 2 p°° _ _ 

Pi = | J H x (y,0)Ey*(y,0)dy = | ^ H x (ky,o)Ey*(k y ,o)dky 


(13) 


( 1*0 


where the dual integral representation is a statement of Parseval's theorem. 
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Using equations (13) and (l4), the external aperture admittance per unit 
length is 


Y Z = 


2Pi 


- i. r bi sln2 ( k y 1) 


e(0) , 

f" g' 

(ky , 0 

) sin2(k y 1) 

e 0 J 

0 ^ 

ky,0 j 

1 1 

^ 1 ) 

1 2 


(15) 


In order to compute radiation patterns, H x ^-(y, z), hence Ej(ky), must he 
evaluated at points far from the aperture. From equations (10) and (l2) , 


M k y) 

Therefore, with the substitutions 


_ y oe ^H +k * Z °) sin(*y |) 
Ll 2 (ky) 


(16) 


v y ? 



z = p cos 0 


y - p sin 0 

H X n (Py0) 

becomes 

Vq Oco 

e “j(ky |+k z z G ^ 

2* J^, 

L 12( k y) 


(17) 


jp(k z cos 0 + ky sin 0) 


dky (18) 


where it has been assumed that 


sin 




= 1 


(19) 


The integral (l8) may be solved by the method of stationary phase. (See 
ref. 1 for application to this geometry.) The result is 
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I 


J ( z o cos 0 + 1 sln 0 ) w v J e (o) k Q cos 0 
= e /V 0 Y 0 <-— - 7v“ 


H x n (p,0) = 




,ry • m 

g hr - sin 0, 


e" Jk ° P (20) 


where the magnitude of the terms in parentheses is the pattern expression, i.e.. 


P(0) = k 0 



( 21 ) 


II. Infinite Slot With TM X Polarization 

The geometry and appropriate field representation is shown in figure 2. 
If the fields are polarized such that the electric field E(y,z) is directed 
along the x-axis the wave equation for E^ is 


dy 2 


+ = 0 


( 22 ) 


As before, Ex^(y, z) and E^^y, z) are defined through the Fourier 
representation, 

Ex 1 = £ f" I’(k y ,*)e Jk y y dk y 

J 

Use of equations (25) then reduces equation (22) to the following equation 
which must he solved numerically 



subject to the starting conditions. 
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’( k y' z o) - 


f ' (ky, z Q ) - J < 


r 

J k o 

\l 

ro 

H 

k y 

-k n J 

i-fef 


o w 

V 

W 

k o 


> 1 


< 1 


> ( 25 ) 


where 




Jk z z 0 


In matrix form, the transformed field solutions at z = 0 are given by 

E x ( k y > °)\ / f ( k y > °) 


Hi 


■(ky,°) 




Er(k y ) 


jk z z 


z^o 


(26) 




Now, if the electric field at the aperture corresponds to a TEp excitation, 
the following transform pairs are defined at z = 0 


E x I (y>0) = E 0 cos(f) 

i x (k y; o) . . *0 . 




k £ y £ IS. 
2 “ J “ 2 


> (27) 


and the aperture admittance per unit length is 

,2 


Yj = — 2E- „ J*fll«b2 r C ° S ‘( ky ll ■ EfecSi 

M 2 J 0 [j, 2 . ( kyt )2] 2 f (V- 0 ) 


dk- 


(28) 


The radiation fields, as derived from a stationary phase evaluation of E x ^ 
is given by 


E x n (p,0) = E 0 e 


« -jk D p 




e ~ - e -J k z z o 


2rtb cos (k Q k sin 0^ 




tlfL = - sln 0 
V k o 


,0^ it 2 - (kob sin 0) 2 


(29) 
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III. The General Aperture 


The development of a general aperture on a ground plane, shown in fig- 
ure parallels the development given above, except a double Fourier integral 
representation is necessary because the fields depend on al 1 three cartesian 
coordinates. This problem was approached by superimposing TE and TM vector 
potentials. By a straightforward application of the boundary conditions, it 
can be shown that the pertinent differential equations are 

_ -\ 


( 30 ) 


With starting conditions 


dff 

dz 2 


+ k r 


ill - |j2 


f = 0 


d% _ l de(z) dg 
a z 2 e(z) dz dz 


+ V 


e(z) 


g = 0 


f (z 0 ) = 1 


f’(z 0 ) = 

fjkojp 2 - 1 

IPI 

> 1 

1^-ko \]l - P 2 

M 

< 1 

g(z 0 ) = 1 




g'( z o) = j < 

pk D n 2 (z o )^0 2 - 1 

1*1 

> 1 

\jkon 2 (z 0 )l/l - 0 2 

1*1 

< 1 


(31) 


where f3 is defined through the transformation k x = koP cos a and 
ky = ko0 sin a. 

Note that equations (30) and (31) are identical to the equations and 
starting conditions which were used for the infinite slot. It can therefore be 
concluded that when the external plasma problem is solved for the infinite slot 
of TE and TM polarization, it is automatically solved for the finite 
aperture . 

The double Fourier transform of the aperture fields may also be set up 
in matrix form defining a transfer function which relates the aperture fields 
to the free— space fields. This relationship is 
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(32) 

If the aperture field corresponds to a TEqi excitation, then 



(33) 


The transforms of equations (33) are 



And, the aperture admittance is 



w 


10 


t 


* 

Equations (52) and ( 5^) may "be used to express Ej and in terms of V Q ; 

therefore, Hy(k x ,ky,o) may also he expressed in terms of V 0 , i.e.. 


f (j. t 0 \ . SByPo J /k x \ /k 0 g\ € (o) f k y\ ( f < \l 

hyl k x ,JCy,U I o\ liT; \T r l~£ + Ik - / lirrJ/ 

v ' k x 2 + ky 2 \\ ° J \ g / 6 o \ K o/ V^Wf 


(56) 


If equations (5*+) and (55) are converted to polar coordinates and substituted 
into equation (55), "the final admittance expression is: 




oo p 2jt 
0 w 0 


cos 

ff- sin 


2 

sin 

(k 0 a 

l-H- cos 
V 2 

a) 

* 2 - 

(k Q b sin a) 

1 

(k Q a \ 

cos a .j 

I 


X /-cos2af^^\ € . ( . ° - + sin 2 a(-^— p dp da 

) \g / e o \ kof /[ 


(57) 


RADIATION FROM SLOT ANTENNAS ON A CYLINDRICAL GROUND PLANE 


IV. The Infinitely Long Axial Slot on a Cylinder 


The geometry is shown in figure 4. An infinitely long slot is cut into a 
long conducting cylinder with H directed along the axis. The structure is 
coated with a plasma, whose dielectric properties vary radially. It is assumed 
that the fields are independent of z so that H z in regions I and II can be 
described by the following Fourier series in azimuth 


H z I (p,0) = Y~ H mz (p)e“j m 0 

m=-oo 




m=-°° 


(58) 


where the radial function, Hmz(p) satisfied the differential equation 

[Hmz = 0 (59) 


1 de dH mz . v 2 

e(p) 

m 2 1 

) ' e d£ ~ + k ° 

e o 

( k op) 2 


11 



By normalizing with respect to C m , equation (39) can "be numerically inte- 
grated provided the boundary conditions are met at p = b. The continuity 
requirement for H z and E0 gives 


hmz(b) = = H m (2)(k Q b) 


hmz ' ("b) - ir 


®m( p) 


dp C. 


= £Cb)fi.(2)* 


m 


■m 


lp=b 


M 


(4o) 


Equations (39) and (^0) then yield normalized field solutions, h m (a) and 
1^' (a) at the aperture. 

Since the azimuthal component of the electric field in the plasma is 
given by 


“< 3 m 0 


¥( p - 2_ 

m=~oo 


1 f 2- * S Taz( p K -ATB$ (4i) 

«<p) fee d P 


the field transforms at the aperture can immediately be written as 


Hmz( a )\ / hmz( a ) \ 


Em0( a ) 




^m 


The power per unit length radiated by the aperture is 

r^ol 2 „ _ 

F l = 2 , E 0 (a^)Hz I ( a ,0)d0 

-0o/ 2 


00 00 


= & 
2 


m=-oo in =-oo 




(^2) 


00 00 


= na Y1 21 V0 I *( a ) H mz I ( a ) 5 m. m ' = rta Ein0 I * H mz I ( a ) 

m=-oo m T =-oo m=-oo 


(43) 


where the equivalence of the integral over physical space to the sum over mode 
space is nothing more than Parseval's theorem in cylindrical coordinates. 
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JJc, is now assumed that the electric field distribution at the aperture is 


E0(a,0) = 12 




whose transform is 


Em0(a) = 


0o/2 
1 2 
fm0 


at J.0 o/2 " 


7 0 


2rta Vm0 o > 


(45) 


From equations (42) and (45) the transform of H z -*- at the aperture is 

Mr 


TT f \ ( \ Vo 2 y t 1 mz (a) 

H mz (a) = -jue(a)— -r, - rrT 

2ita (m</) 0 \ hj^ (a) 


Therefore, the admittance per unit length is: 

00 

XL 


Y 2P .koYo eu 

1 " lv„l 2 ' ' J at B € 


o m=-°° 


' m 




sin 


Mo\ 
\ 2 J 


bm(a) 

4m* ( a ) 


(**6) 


(47) 


The radiation pattern expression is found by solving for substituting 

into the second of equations ( 58 ), and evaluating H z n (p,0) in the far field 
(i.e., by asymptotically expanding the Hankel function). From equations (^2) 
and (45), 


c -■*«(») 7 
^ 2rtah m '(a) ° ^m0 o ^ 


Therefore, H z ^(p,0) in the far field is given hy 

,- J ( k ° P "4kol 


(48) 




Jtkop 


r 

m=-°o 


The pattern is computed by taking the magnitude of the term in brackets 


L#)l 

.-H* -Dl 

t 

hm‘ ( a ) ( 

ro 

J 

the term 

in brackets 


(49) 
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V. The Arbitrary Aperture on a Cylindrical Ground Plane 

The preferable method of approaching the problem of a general aperture 
radiating into a radially varying plasma as shown in figure 5 would be to con- 
struct the fields from axial TE and dM components and numerically solve the 
differential equations. However, such an attempt fails because the Maxwell 
equations will not be satisfied. 

As an alternative approach (ref. 16), the inhomogeneous plasma was 
replaced by several homogeneous tandem layers judiciously chosen to fit the 
given distribution. The development of the problem follows from a straight- 
forward generalization of Wait (ref. 3) , and will therefore only be outlined 
here. 


Since each layer is homogeneous, the axial components of E and H com- 
pletely specify the external problem, and can be written as follows: 


E z 1 (p,0,z) = e" J u i 2 F m 1 (k z ,p)e _JkzZ dE z 


i m=-o° 


} (50) 


00 / O 00 

H z 1 (p,0,z) = — i S'" e -jm0 / u i 2 G a i (k z ,p)e“'5 k 2 z dk z 

frVoH m=~ J X 


where the superscript i refers to the fields in the i^* 1 layer. When 
i = 1, 2, . . ., n, the radial solutions of equations ( 50 ) are: 


Em 1 (k z ,p) = a m 1 J m (u i p) + A^u-^) 

G m 1 (k z ,p) = ’ b m ij m (uiP^ + B^Uip)! 


(51) 


Wo standing waves exist in free space (i = n + l) ; therefore, F m ^ and Gm 1 
are of the form: 


F m n+1 (k z ,p) = c m H m ( 2 )(u 0 p) 

Om m+ ^ (^z > p) = ^mAn ^ ^ p) 


(52) 


If equations (50) are used in connection with the Maxwell equations, the 
azimuthal components of E and H can be computed, and the boundary condi- 
tions can therefore be stated at each interface. From the schematic shown in 
figure 6, the boundary conditions at p = P±+± are: 


Ik 



E 0 1 ( p i+l , ^ ,z ) “ E 0 1+1 ( p i+l^ ,Z ) 

E z 1 ( p i+l>^ z ) 

H 0 1 ( p i+l^ ,z ) 

H z ^ (pi+ 1^0^ z) = H z 1 + 1 (p^ + -p 0 ,z^j 

The fields are defined in terms of the unknown coefficients, a^, b m , Am, and 
Bjq; therefore, equations (53) can he expressed in terms of 4 X 4 matrices 
relating the unknowns in the i 1 ^ 1 layer to those in the (i + l) 1 * 1 layer, i.e.. 





where the A and B matrices are given explicitly in appendix III. 


If another 4x4 matrix is defined such that 



then it can easily he verified that 



That is, the unknown coefficients in the first layer can he related directly to 
those in the last layer hy multiplying similar 4x4 matrices. 


If the boundary conditions are completed at the air -plasma interface, and 
if am 1 , Am 1 , bm 1 , and Bm 1 are converted to field transforms at p = a, the 
following transfer function is defined 
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(56) 


r -> 

Emz ^a,k z ^ 


r ~\ 

0 

H m 0^a,k z j 
Hmz( a ^^z) 

) ~ |^jk( m > k zj] 

Cm^z) 

0 

J^n0( a > k z) 

1 

^mW 

V ) 


where the 4 X 4 L-matrix is entirely composed of specified or computable param- 
eters (dielectric constant, Bessel functions, etc.). The L-matrix is symboli- 
cally defined in Appendix III. The admittance and patterns can be found by 
completing the boundary conditions at the aperture, and algebraically solving 
for the unknowns. 


RESULTS 


In order to compare the radiating characteristics of slots on flat and 
cylindrical ground planes under typical reentry conditions, the electron den- 
sity and collision frequency profile shown in figure 7 was chosen as a coating. 
The slot on both ground planes was infinitely long, 0.4 inch wide (the -width 
of a standard X-band waveguide), and polarized with the magnetic field directed 
along the slot. The patterns and input admittance were then computed via the 
methods of sections I and TV. 

The radiation patterns of the coated slot on a flat ground plane are shown 
in figure 8 . The behavior of the patterns as a function of frequency is sharply 
divided at the peak plasma frequency of 10 . 76 gc. Frequencies above this 
exhibit definite off-axis maxima which progress toward the axis as the frequency 
decreases. Whereas, below 10 . 76 gc, the maxima occur only at 0=0, and the 
radiation rapidly decays as a function of angle. The shifting of the maximum 
and the differences in the pattern characteristics on either side of critical 
frequency are exactly what one would expect if the coating were homogeneous with 
a plasma frequency of 10 . 76 gc. However, the location of the pattern maximum, 
which for a collisionless homogeneous plasma is given by the plane wave critical 
angle 


0 C = cos - 1 (fp/f ) 

does not correlate with the results of figure 8 . 

The radiation patterns of the coated cylinder are shown in figure 9* 
Again, the radiation patterns have strikingly different characteristics, 
depending upon whether the propagating frequency is greater than or less than 
the peak plasma frequency. There is also a progression of the pattern maxima 
above plasma resonance, but the location cannot be predicted as easily as the 
plane case. 
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• * The external admittance of the cylinder is plotted as a function of fre- 
quency in figure 10. Nineteen admittance points were computed in the frequency 
interval of 10.0 to 11.8 gc, yet this number was insufficient to establish a 
smooth curve because of small-scale fluctuations; nonetheless, some interesting 
features are revealed. The most striking effect occurs in the region of peak 
critical frequency, where the conductance and susceptance suddenly decrease. 
Based upon some recent work by C. M. Knop, low values of admittance should be 
expected at critical frequency. 

At frequencies above resonance, the conductance rises rapidly, but the 
susceptance remains relatively constant. As the frequency increases, both 
curves smoothly approach the no -coating values, as they should. Below reso- 
nance the curves are fluctuating too much to make any general comments. 

In order to save computer time, the plane results were confined to the 
five points shown in figure 10. The computations are not as extensive as those 
for the cylinder; nevertheless, there are sufficient data to conclude that for 
the parameters considered here, the admittance of a slot on a plane differs 
from the cylinder by only a few percent over the entire frequency range. 

If the losses are small (as they are for this plasma distribution), trapped 
waves in azimuth are supported by the cylindrical structure at frequencies above 
resonance. Therefore, it is conceivable that the wave interference within the 
cylindrical coating could be such that the difference between the admittance of 
this structure and that of a slot on a flat ground plane would be substantial. 
The results of figure 10 seem to indicate that this is not so, at least for the 
size cylinder considered here. It is of interest to note, however, that the 
greatest differences in the admittance occur at frequencies above resonance, 
which confirms the intuitive expectation. 

The computations of the external admittance of a finite aperture, coated 
with the inhomogeneous plasma, have not been done; however, the ratios 


s (Ml 

g'(P,o) 


and 


f'(p,o) 

f(p,0) 


appearing in the integrand of equation (37) are given as a function of 3 in 
tables 1(a) through 1(d) . Using these tables, the double integral which defines 
the admittance can be numerically evaluated. 

The equatorial (0 = jt/ 2) radiation patterns of a half-wave axial slot, 
coated with a homogeneous plasma, are shown in figure 11, and were derived by 
two independent methods. The points were derived from the bracketed term of 
equation (^9), 2 and the solid curve was derived from the L-matrix of section V. 
The explicit expression for equatorial pattern, in terms of the L-matrices from 
the stationary phase evaluation (ref. 7) is 


^In the equatorial plane, the patterns of a finite slot are identical to 
those of an infinitely long slot (ref. 7) • 
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’(H - 




Tt c 


00 

\y 

j m cos m0 

z 

1 m=0 

(1 + 


The agreement is good; therefore, one method substantiates the other. 

The elevation patterns shown in figure 12 were taken directly from the 
results of W. V. T. Rusch (ref. 6), who based his treatment on a thin sheath 
approximation. The use of this approximation assures that the solution within 
the plasma may be expanded in a Taylor series about the thickness b-a. How- 
ever, as the polar angle decreases from 90 ° to 0°, this approximation becomes 
invalid. The points on figure 12 were computed from the expression^ 


( Q o\ _ 

2v3Lq cos^i cos 0^ 


OO 

y 

« 5 ±n ‘ L 12 

V ’«) ~ 

rt^sin 0 


z 

m=0 

(l + 5 0 m ) (L 12 Ll^ - L4 2 Li4] 


and were normalized with respect to Rusch 1 s on-axis value of mjay p = 1.048 
curve. There is reasonably good agreement at 0 = 70° and 0 = 90°, "but the 
results begin to deviate as 0 decreases. However, this was expected, so it 
is reasonable to assume that the method can validly be extended to the admit- 
tance problem. 


CONCLUDING REMARKS 


Methods have been developed whereby the radiation patterns and near-field 
effects can be evaluated by solving the aperture and plasma problems inde- 
pendently. As a special application of these techniques, infinitely long slots 
on planes and cylinders, coated with a reentry plasma were analyzed, and the 
following conclusions were noted: 

(1) If the collision frequency is low, the peak value of the plasma fre- 
quency determines the radiating characteristics of the aperture. At critical, 
both the patterns and admittance undergo drastic changes. 

(2) A homogeneous plasma " equivalent" of the inhomogeneous plasma coating 
given in figure 7 will give only general trends of pattern behavior. 

( 3 ) The external admittance of the particular slotted planes and cylinders 
considered here differs by only a few percent. 

Tables were also given, whereby the external admittance of a finite aper- 
ture on a plane coated with the plasma of figure 7 can be computed by using 
these tables in connection with the double integral (37). 


5rhe positive or negative sign on j is chosen if 0 is 0 or m, 
respectively. 
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•» • The solution finite aperture on a cylinder was formulated in terms of 
matrix products and verified hy comparing the radiation patterns with known 
results. 
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APPENDIX I 


LIST OF SYMBOLS 


a 

b 

E 

E,E 

Eo 

f 

f P 

H 

H,H 

Ho 


slot dimension along the x-axis or radius of cylinder, as appropriate 

slot dimension along the y-axis or radial distance to air-plasma 
interface, as appropriate 

electric field intensity 

Fourier transforms of the electric field 

electric field intensity at aperture of TM X slot 

propagating frequency 

plasma frequency, 8.97 X 105 Jn^" 

magnetic field intensity 

Fourier transforms of magnetic field 

field intensity of incident TEM mode in waveguide 


j = fi 


k Q free-spaee wave number 

k x ,ky,k z Cartesian mode numbers 


m 

n 


N e 

P 

r,e,0 

u 

Vo 

x,y,z 


azimuthal mode numbers 

index of refraction or discrete mode number in waveguide, as 
appropriate 

electron density- 

power 

spherical coordinates 
radial mode numbers 
applied potential on slot 
Cartesian coordinates 


20 


admittance of free -space, 


Y 

Zo 


admittance 


impedance of free-space. 



z Q plasma thickness 

ky 

a = arc tan — 

k x 

P = l|k x 2 + k y 2 

r reflection coefficient 

8 m n Kronecher delta (m = nj 

e permittivity 

e Q permittivity of free space 

Iq free-space wavelength 

permeability of free space 
v collision frequency 

p, 0 , z cylindrical coordinates 

0 C plane-wave critical angle 

0 q angular width of slot 


0 angular frequency 
* complex conjugate 
Superscripts : 

1 fields within the plasma 


n 


fields in free space 


i fields within the i^h dielectric layer * ' 

n fields within the n^h dielectric layer 

Subscripts: 

i physical parameters within the i th dielectric layer 

l per unit length 

x,y,z vector components in the principal Cartesian directions 

p,0,z vector components in the principal cylindrical directions 

A prime denotes differentiation with respect to z or p, as appropriate. 
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APPENDIX II 

EXACT FORMULATION OF THE INFINITE SLOT ON THE GROUND PLANE 
RADIATING INTO AN INHOMOGENEOUS PLASMA 


Assuming a TEM wave incident upon the aperture shown in figure 1, the 
pertinent internal fields at z = 0 are given by 


H x (y,0) = Ho 


1 + XI A n cos (ir) 

n=0 ' D L 


-stes'-S [>-£«■ ™(T) 


(II-l) 


where 




- fif - Wf 

■ 3 " ''o 2 

The external fields at the aperture are, from equation (10) 

H x(y»°) - sj / Lai(iv) B T e3kzZoeJ%yak y 

^ —OO 

/ oo 

i a2 (v) B r eJkzl0eJk3rydi y 


hu >m 


mr 


> (II-2) 


kn < 


J 


(II-3) 




The continuity requirement on Ey gives 

k f £ ( 8 ° n - *$£ AT ^ 


Multiplying equation (II-4) by e“*% y dy, and integrating with respect to y 
gives 
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(II-5) 


L 12 (ky)nTe Jk ^° = 1 — 


o n=0 



COS 2S e - Jk y y dy 

b 


where the integration of the left-hand side of the above was extended over all 
space because Ey = 0 on the ground plane. 

It is now convenient to define 


1 

b 



cos 221 e +;|k y y ay ~ jk y b .. - R-1)V¥ - ll 

b (k y b)2 . ( njI )2 L J 



(II-6) 


which has the property 



c n( k y) c Z*( k y)<% = fn l (l + & 0 n ) 


(II-7) 


that is, C n ^ky^ is an orthogonal function. Therefore, in terms of C n (ky), 
the continuity of E gives the following equation 


e J k z z o 

Ho 




(II-8) 


The boundary conditions for H x is a little more difficult to apply 
because the magnetic field involves "mixed" boundary conditions, i.e., H x is 

continuous at the aperture, and equal to the unknown surface current on the 
ground plane. The boundary condition must therefore be applied only over the 

region where the aperture exists. This is done by expanding e^ k y y in a 
cosine series over the aperture; i.e., 

e Jk y y = B n cos 22X 0 5: y 5? b (ll~9) 

— r\ b 


Use of orthogonality gives 


B n 



(II-10) 


Therefore, 
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(11-11) 


e j%y _ g , c n (k y ) 

fco (1 + s 0 n ) nV y ' 


COS 


nity 


Hr 


or 


Use of equation (il-ll) in the boundary condition for H x , gives 

oo poo oo 

5 (V + A n )cos H.l/^ ^ 

( 8 °" + ^ = n (l + S 0 n ) £ eJkz2 ° C n(%) a V 


COS 


nity 


b 

( 11 - 12 ) 

(11-15) 


By substituting equation (II-8) into equation (II-13), 


(b Q n + A n ^ = 7— ^ -fpr f 

V ' Y 0 * (l + B 0 n ) J -c 


dk-u 


L n( k y; 


C n(*y) £. - A q ): 


qz n * 


( k y) 


(n-14) 


Now, let 


J nq “ 


r(l + »o n ) 

The equation (11-1*0 becomes: 


’ Y 0 ,(l\ So n\ ^ (II ‘ 15) 


V CO 

K + A n ) = ^ (8o q - A,) 


nq 


or, 


(11-16) 


) (inq + = ^no “ ^n 

q=0 ' ' 


(11-17) 


Therefore, since the coefficient matrix, J n ^ can be computed after L]j_ and 
Lip are specified, the problem reduces to solving an infinite matrix (II-16) 
in order to compute the unknown A^'s. 
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APPENDIX III* 4 - 

THE FINITE APERTURE RADIATING INTO 
LAYERED DIELECTRICS 


Upon applying the boundary conditions (55) at p = p i+ ^_, the matrix 

relationship between the coefficients in the i^* 1 and the (i + l) s t layer can 
be expressed in many ways; however, if closed form expressions are desired, 
the most profitable way seems to be 


[V +1 ] 


J m(uiPi+i) 

u iV ( u iPi+l) 

0 


hi 


(^ 1 %) j “( u i+i p i+i) 
“l+y®' ( u i+i c i+i) 


H m (2) (u iPi+1 ) 

u i H m^ ' ( u iPi+l) 


“i+lV 2 ’ ' (“l+lOl+l) 


0 


I ^Cpink^ 2 / kj 2 

k o 2 u i 2 Pi+l\ *i+l 2 ) 






( u ipi+:) 


u i J m' ( u iPi+l) 


0 


H m (2)( u i p i+1 ) 

’ ( u iPi+l) 


(HI-1) 


>up 0 mk z ^ _ k 
u i 2 °i + i 


: i 

k i+l‘ 


^ u i+lPi+lJ 


i j i.k + ipi + i) - i^?) B - <2) ( u i + i p i + i) 

tetf Jn(ui+iPi+i ) te£f h " <2) 

V K+l D l + l) Ul+1 ( E ltiJ °”* 2) ^ 1+1 ° 1+1 ) 


(III - 2 ) 


where the 


A jk i+1 


L 


and 


C matrix was defined as 


V +1 ] 

Vi 


matrices are defined in equation (5*0 • The 



(HI-5) 


^Excerpts from reference 16 


were used in this appendix. 
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where the inverse of the A-matrix is given by 


[ V * 1 ]' 1 - 


j| u iPl+l H m (2) Kp 1 + x) -j|pi + l H m (2) Kp 1 + i) 
"•52 U l p i+l Jm ( u l p i+l) J? p i+l J m( u iPi+l) 


0 0 (^Pi+l) -J | p i+l H ni^ 2 H u i p i+l) 

0 0 -Jf^Pi+lJm' ("lPi+^ jfp i+1 Jm(uiP i+ l) 

(III -4) 

If equation (55) is carried one step further, it follows that 


i) 

a m 


r -\ 
0 

Am 1 

>■ - [ c Jk 2 ][ c Jk 5 ] • • • [ c jk n - 1 J[fj k n )[cjk n+1 ] < 

0 

yj 


J 


where use was made of the fact that 


C +1 ■ o 


A m n+1 = % 


V + 1 - 0 


B m n+1 - 3m 




(III-5) 


(III-6) 


describe outgoing waves in free space. From the explicit form of the field 
transforms, it can be shown that 
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(Ill -8) 


Hence, the L-matrix of equation ( 56 ) is the product of the n + 1 matrices 


L jk - ^k 1 c jk 2 c jk^ 


r n P n+1 
C jk c jk 


It is interesting to note that if the index of refraction of all the layers is 
unity, the C-matrices are unit matrices. 
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g(P,o) 

g'(P,o) 


f'(P,o) 

f(P,o) 


p 


0 

0.035 + jo.163 

- 0.424 - j 0.844 

0.050 

0.034 + j 0 .i 64 

- 0.427 - jo .842 

0.100 

0.030 + jo.166 

-0.433 - jo.836 

0.150 

0.023 + jo. 168 

-0.444 - jo.825 

0.200 

o.oi 4 + jo.171 

- 0.460 - j 0 . 8 u 

0.225 

0.009 + j.0172 

- 0.469 - jO.803 

0.250 

0.003 + jO.174 

- 0.480 - jO -793 

0.275 

- 0.004 + jo.175 

-0.492 - jo.783 

0.300 

-0.011 + jo.176 

-0.504 - jo.771 

0.325 

-0.018 + jo.176 

-0.518 - j 0.759 

0.350 

-0.026 + jo.176 

-0.55^ - , 30.746 

0.375 

-0.034 + jo.176 

-0.550 - jO.752 

0.400 

- 0.042 + jO.175 

-0.568 - jO.716 

0.425 

-0.050 + jo.173 

-0.586 - jO.700 

0.450 

-0.058 + jo.171 

-O.606 - j 0.684 

0.475 

-0.067 + jo.169 

-0.627 - jO.666 

0.500 

-0.074 + jo.166 

-0.650 - jo .647 

0.550 

-0.089 + jo.159 

-0.698 - jO.608 

0.600 

-0.103 + jO.150 

-0.752 - jO.565 

0.650 

-o.n 4 + jo.i 4 i 

-0.810 - jO.520 

0.700 

- 0.124 + jo.131 

-0.874 - jO.471 

0.800 

-0.138 + jO.lll 

- 1.016 - jO .364 

0.900 

-0.147 + jo.093 

- 1.180 - J 0.240 

1.000 

-0.152 + jo.082 

- 1.364 - J 0.002 

1.200 

-0.156 + j 0.048 

-2.027 - jO.OOl 

1.400 

-0.150 + jO.031 

-2.508 - jO.OOl 

1.600 

- 0.144 + j 0.021 

-2.966 - jO.OOl 

1.800 

-0.138 + jO.Ol 4 

- 5.413 - jO.OOl 

2.000 

-0.133 + jO .010 

-3.855 - jO.OOl 

2.500 

- 0.121 + j0.005 

- 4.946 - jO.OOO 

3.000 

- 0.110 + j 0.002 

- 6.028 - jO.OOO 

3.500 

- 0.101 + jO.OOl 

-7.102 - jO.OOO 

4.000 

-0.092 + jO.OOl 

-8.171 - jO.OOO 













p 


g(P,o) 

g'(P,0) 


f*(P,o) 

f(p, 0 ) 


0 

0.070 + JO. 200 

-0.280 - JO. 985 

O.lOO 

0.074 + JO. 250 

-0.290 - JO. 975 

0.200 

0.055 + j 0.549 

-0.520 - JO. 945 

0.250 

- 0.020 + JO . 446 

- 0.542 - JO. 921 

0.500 

-0.188 + JO. 498 

-0.570 - JO. 894 

0.550 

- 0.564 + J 0.410 

-0.405 - JO. 865 

o. 4 oo 

- 0.422 + jO.261 

- 0.442 - JO. 827 

0.450 

- 0.402 + JO. 155 

- 0.486 - JO. 787 

0.500 

-0.565 + JO. 095 

-0.556 - JO. 745 

0.550 

-0.544 + JO. 070 

- 0.469 - JO. 715 

0.570 

-0.520 + JO. 050 

-0.615 - JO. 676 

0.600 

- 0.504 + JO. 059 

-0.655 - jo. 645 

0.650 

-0.291 + JO. 052 

-0.695 - JO. 612 

0.670 

- 0.275 + j 0.024 

-0.749 - JO. 567 

0.700 

- 0.266 + JO . 020 

-0.794 - JO. 552 

0.725 

-0.259 + jo. 017 

-0.855 - jo. 502 

0.750 

-0.252 + jo. 015 

-0.875 - jo. 471 

0.775 

- 0.246 + jo. 015 

-0.916 - J 0.440 

0.800 

- 0.241 + JO. Oil 

-0.959 - JO . 407 

0.820 

-0.257 + JO. 010 

-0.995 - JO. 580 

o. 84 o 

-0.255 + jo. 009 

-1.052 - JO. 555 

0.860 

-0.250 + JO. 008 

-1.070 - JO . 524 

0.880 

-0.226 + JO. 008 

-1.109 - JO. 295 

0.900 

-0.225 + jo. 007 

- 1.149 - jo . 264 

0.920 

- 0.221 + JO . 006 

-1.190 - JO. 252 

0.940 

-0.218 + JO. 006 

-1.252 - jo. 197 

0.960 

-0.215 + jo. 005 

-1.275 - JO. 157 

0.970 

- 0.214 + JO . 005 

-1.297 - JO. 155 

0.980 

-0.215 + jo. 005 

-1.519 - JO. 109 

0.990 

- 0.212 + JO . 005 

-1.541 - JO. 077 

0.995 

- 0.211 + J 0.004 

-1.552 - JO. 055 

1.000 

- 0.210 + J 0.004 

- 1.564 - JO . 002 

1.200 

- 0.191 + JO . 002 

- 2.102 - JO . 001 

1.500 

-0.172 + JO. 001 

-2.895 - JO. 001 

2.000 

- 0.148 + JO . 001 

- 4.121 - JO . 000 

2.500 

-0.129 + JO. 000 

-5.515 - jo. 000 

5.000 

-o.n 4 + jo. 000 

- 6.488 - J 0.000 

5.500 

- 0.102 + JO . 000 

-7.654 - -JO. 000 

4.000 

-0.092 + JO. 000 

-8.812 - JO . 000 







TABLE I.- TABLES OF 


g(P,Q) 

g'(P,0) 


AND 


FOR THE INHOMOGENEOUS 


f(0,O) 


PLASMA GIVEN IN FIGURE 7 - Concluded 
(d) f = 11.4 gc 


0 

g(P,0) 

g'(0,O) 

f'(0,O) 

f(P,0) 

0 

0.050 + 30.194 

-0.168 - 31.099 

0.050 

0.051 + 30.196 

-0.170 - 31.096 

O.lOO 

0.055 + J0.202 

-0.178 - 31.087 

0.150 

0.057 + J0.212 

-0.192 - 31-075 

0.200 

0.041 + 30.228 

-0.211 - 31.055 

0.225 

0.044 + 30.258 

-0.222 - 31.041 

0.250 

0.046 + 30.251 

- 0.255 - 31.028 

0.275 

0.048 + 30.265 

-0.250 - 31-015 

O.JOO 

0.050 + 30.285 

-0.266 - 30.997 

0.525 

0.050 + 30.504 

-0.285 - 30.979 

0.550 

0.049 + 30.529 

-0.502 - 30.960 

0.575 

0.044 + 30.558 

-0.522 - 30.940 

o.4oo 

0.055 + 30.592 

-0.544 - 30.919 

0.425 

0.019 + 30.451 

-0.567 - 30.897 

0.450 

-0.008 + 30.474 

-0.592 - 30.875 

0.475 

- 0.049 + 30.518 

-0.4i8 - 30.849 

0.500 

-0.108 + 30.558 

- 0.446 - 30.825 

0.550 

-0.277 + 30.588 

-0.507 - 30.768 

0.600 

-0.449 + 30.502 

- 0.575 - 30.710 

0.650 

-0.529 + 30.546 

- 0.649 - 30.648 

0.700 

-0.520 + 30.210 

-O.750 - 30.582 

0.800 

-0.429 + 30.071 

- 0.914 - 30.441 

0.900 

-0.555 + 30.024 

-1.125 - 30.285 

1.000 

-0.505 + 30.002 

-1.564 - 30.002 

1.200 

-0.259 + 30.001 

-2.160 - 30.001 

l.4oo 

-0.208 + 30.001 

-2.757 - 30.001 

1.600 

-0.187 + 30.000 

-5.279 - 30.001 

1.800 

-0.170 + 30.000 

-5.805 - 30.000 

2.000 

-0.157 + 30.000 

-4.522 - 30.000 

2.500 

-0.155 + 30.000 

-5.589 - 30.000 

5.000 

-0.115 + 30.000 

-6.856 - 30.000 

5.500 

-0.101 + 30.000 

-8.069 - 30.000 

4.ooo 

-0.091 + 30.000 

-9.294 - 30.000 













Figure 3.- Aperture coated with inhomogeneous plasma. 





DIELECTRIC 



Figure 5*- General aperture on a coated cylinder. 







Figure 6.- Schematic for applying boundary conditions. 




Figure 8.- Patterns of the coated slotted plane. 




Figure 9*- Patterns of the coated slotted cylinder. 




Figure 10.- Admittance of coated slots. 



Figure 11.- Equatorial patterns of coated axially slotted cylinders 






